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NOTE: This paper presented the first attempt to tackle the Kardar-Parisi- Zhang (KPZ) equation 
using non-perturbative renormalisation group (NPRG) methods. It exploited the most natural and 
frequently used approximation scheme within the NPRG framework, namely the derivative expansion 
(DE). However, the latter approximation turned out to yield unphysical critical exponents in dimen- 
sions d > 2 and, furthermore, hinted at very poor convergence properties of the DE. The author has 
since realized that in fact, this approximation may not be valid for the KPZ problem, because of the 
very nature of the KPZ interaction, which is not potential but derivative. The probable failure of the 
DE is a very unusual - and instructive - feature within the NPRG framework. As such, the original 
work, unpublished, is left available on the arXiv and can be found below. 

Added note: the key to deal with the KPZ problem using NPRG lies in not truncating the momen- 
tum dependence of the correlation functions, which is investigated in a recent work arXiv:0905.1025 . 

We present a new approach to the Kardar-Parisi- Zhang (KPZ) equation based on the non- 
perturbative renormalisation group (NPRG). The NPRG flow equations derived here, at the lowest 
order of the derivative expansion, provide a stable strong-coupling fixed point in all dimensions d, 
embedding in particular the exact results in d — and d = 1. However, it yields at this order 
unreliable dynamical and roughness exponents z and \ in higher dimensions, which suggests that a 
richer approximation is needed to investigate the property of the rough phase in d > 2. 

PACS numbers: 64.60.Ht,05.10.Cc,68.35.Fx,05.70.Ln 

The Kardar-Parisi-Zhang (KPZ) equation, though originally introduced as a coarse-grained de- 
scription of non-equilibrium interface growth [l[ , has acquired a broader significance over the past 
decades as a simple model for generic scale invariance and non-equilibrium phase transitions 0. 
It is indeed intimately related to many other important physical systems, such as randomly stirred 
fluid (Burgers equation) 0] , directed polymers in random media [1] , dissipative transport [(| 
or magnetic flux lines in superconductors 0]. The KPZ equation has thus emerged as one of 
the fundamental theoretical models to investigate universality classes in non-equilibrium scaling 
phenomena and phase transitions [2[ . 

It is a non-linear Langevin equation which describes the large- distance, long-time dynamics of 
the growth process specified by a single- valued height function /i(x, t) on a d-dimensional substrate 
x e M d : 

d t h(x,t) =j/V 2 /i(x,i) + A/2 (V/i(x,t)) 2 + T){x,t), (1) 

where r/(x, t) is a zero mean uncorrelated noise with variance (f?(x, t)rj(x' , t')) — 2£><5 d (x — x') S(t — 
t'). This equation reflects the competition between the surface tension smoothing force v\/ 2 h, the 
preferential growth along the local normal to the surface represented by the non-linear term and 
the Langevin noise r\ which tends to roughen the interface and mimics the stochastic nature of the 
growth. 

The objective is to determine the profile of the stationary interface, characterised by the two- 
point correlation function C(|x— x'|, t — t') = ([h(x., t) — h(x' , t')] 2 } and, in particular, its large-scale 
properties where C is expected to assume the scaling form C(L,t) = L? x f(r/L z ), where \ and z 
are the roughness and dynamical exponents respectively. These two exponents are not independent 
since the Galilean symmetry Q — the invariance of Eq. (JXJ) under an infinitesimal tilting of the 
interface — enforces the scaling relation z + x = 2. 

While the KPZ equation is solvable in d = 1 [J, Q and on a Bethe lattice 0], its complete 
theoretical understanding is still lacking. The one-dimensional case is special due to the incidental 
existence of a fluctuation-dissipation theorem (FDT) which yields the exact result \ — 1/2 and 
z = 3/2. For d > 2, Eq. JT]) entails two different regimes, separated by a critical value A c of 
the non- linear coefficient [U, y]. In the weak-coupling regime (A < A c ), the behaviour is governed 
by the A = fixed point — corresponding to the linear Edwards- Wilkinson equation 0] — with 
exponents \ = (2 — d)/d and z = 2. In the strong-coupling regime (A > A c ), the non-linearity 
becomes relevant and despite considerable efforts, the statistical properties of the rough interface 
above d = 2 remain very controversial. 
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Indeed, few analytical approaches have proved successful. While the numerical integration of 
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Id . perturbative RG, even re-summed 
On the other hand, surprisingly 



discretised versions of Eq. (JTJ) is plagued with singularities 
to all orders, fails to describe the strong-coupling phase 

accurate exponents can be worked out from mode-coupling theory [l2j but the latter would 
deserve to be put on firmer grounds since it remains a 'ad- hoc' approximation ||. As a consequence, 
the very existence of an upper critical dimension d c for this problem is still much debated. Most 
of the analytical arguments [ill, El, Ell suggest d c ~ 4, whereas numerical simulations and real 
space calculations find no hints for a finite d c [14| . 

In this context, a controlled analytical approach appears highly desirable. This motivates our 
resort to the NPRG method [H, [l6[ . This formalism, which has allowed to tackle notoriously dif- 
ficult str ong -coupling problems at equilibrium [l5| , has been recently extended to non-equilibrium 
systems [IlCiJ]. It has given rise to important progress for reaction-diffusion processes, revealing 
crucial non-perturbative effects [13, El] • It therefore stands as a promising candidate to investigate 
the KPZ problem and we show in this work that it indeed provides the strong-coupling fixed point 
that allows to describe the rough phase in d > 2. 

Our starting point is the field theory corresponding to Eq. |T]), with response functional: 



S[h,h] = J d d xdt Mi 



d t h-vV 2 h-^ (Vh) 2 



Dh 2 \, (2) 



which follows from the Janssen - De Dominicis formalism [l9[ . This field theory can be investigated 
using the NPRG. We don't detail here its implementation but rather set out its principle (see [l5| 
for reviews). This formalism relies on Wilson's RG idea, which consists in building a sequence of 
scale-dependent effective models, that interpolate smoothly between the short-distance physics at 
the (microscopic) scale k — A and the long-distance one at k = 0, through progressively averaging 
over fluctuations. Rather than expressing — as in the original Wilsonian formulation — the flow 
of effective Hamiltonians for the long-distance modes, one can work out the flow of effective 'free 
energies' V k for the short-distance modes that have been integrated, following [l5j]. Thus, at 
k = A, no fluctuation has yet been included and Ta coincides with the microscopic action S, while 
at k = 0, all fluctuations have been integrated out and Tq is the analogue of the Gibbs free energy 
r at thermal equilibrium, in that it encompasses the large-scale (and long-time) properties of the 
system. To construct Tfc, one needs to suppress the low-momentum modes. This is achieved by 
adding a scale-dependent term to the original action [la, |2Cj : 

AS k [h,h]= [ q 2 M-q,- W )i?-(q 2 )fo(q, W )-2/ l (-q,- W )i?°(q 2 )/ l (q,c) (3) 

where the cutoff function i?£ (q 2 ) = X k r(q 2 /k 2 ) is large for small momenta |q| < k so that the 
low-momentum modes are decoupled and vanishes for large momenta |q| > k. The generating 
functionals Zk[j,j] = j T>hT>ih exp(— S — ASk + J jh + J jh) become therefore fc-dependent. Tfc is 
the (modified) Legendre transform of log Zk\j,j]: Tk[ip,tp] = — log Z k [j, j] + J jip+ £j]Jj—ASf.[ip,ip], 
where the last term ensures that Tfc has the proper limit at k = A: ~ S [l5l . Hc|. Tj, is a 

functional of the conjugate fields ip = 8logZk/Sj and ip = S\ogZk/Sj. An exact functional 
differential equation governs its RG flow under an infinitesimal change of the scale s — log(fe/A) 
51123: 



;)J,=~d a Tr[ lnP[V>,# (4) 



where d s = d s Rf! d^D + d s R^ and P[ip, ip] is the 2x2 matrix of second derivatives of (r& 
+ ASk) with respect to (w.r.t.) ip an d ?P- P 1 embodies the full ( i.e. functional) renormalised 
propagator of the theory. Obviously, Eq. cannot be solved exactly and one usually truncates 
it. However, as the truncations used do not rely on the smallness of a parameter, the approach 
remains in essence non-perturbative. 

Since the critical physics corresponds to the lon g-di stance, long-time limit, a sensible truncation 
consists in expanding Tk in powers of gradients [la ] and time derivatives. At leading order in 
derivatives, the Ansatz for Tk related to the field theory ^ writes: 

T k = J ' d d xdt{i>[ti k d t -v k V 2 ]il;-~\ k >4)(Vi)) 2 -D k $ 2 . (5) 
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This Ansatz must be invariant under a Galilean transformation, which writes: 



r ^ 



V»(x, t) — ► V( x + Afcvt,t) + v.x 



(6) 



Moreover, the special time reversal symmetry that yields FDT in d = 1 can be encoded as an 
invariance under the field transformation 



r -0 



ip(x-,t) -> -V>(x, -i) 
i^(x,t) - ^( x ,-t) + ^-A^(-t) 



(7) 



which is satisfied both by the bare action (J2J) and by our Ansatz (J3J). 

First, note that A5fc[^, V 7 ] defined by (|3j has been deviced to respect both the Galilean and time 
reversal invariance. Then, the requirement that ([5]) is Galilean-symmetric implies fj, k A& = Xk , that 
is /ik = 1) which means that fields are not renormalised — they acquire no anomalous dimension. 
This feature can be checked explicitly by calculating the flow equation for fi k — one finds d s pk = 0. 
Also, Afe, structure constant of the Galilean symmetry group, should not be renormalised (as can 
be inferred from the corresponding Ward identities |21|). We show below that this is indeed 
automatically satisfied in our formalism. 

The non-perturbative flow equations associated with the Ansatz ([5]) can be derived as follows. 
The coupling constants Vk, Xk and Dk are defined from ((SJ) by: 



Vk = 9 a 



d 2 T k 
dip dip ' 



A fc = d„2 ^ 



1 d 2 T k 



$ dip 2 



and Di. = — — 



ld 2 T k 



2 dip 2 



(8) 



evaluated at a uniform and stationary field configuration {ipo 1 tpo). Their flow equations follow from 
Eq. Q by taking the corresponding derivatives w.r.t. fields (and momentum), and evaluating them 
at (ipo,ipo). After integrating over the frequency u>, we obtain: 



d s D k 



\ 2 



^{Dk + R D k \ 



(2tt)° 



q 4 (v k + R v k ) 2 - 2 q 2 X k 4>o {D k + R%) 



3/2 • 



(9) 



d a Vk 



A 2 , d d q 4 (2-d)(C fe + i?-)K + i?-) + 2 q 6 ((C fe + ^)a q2 ^-K+^)9 q2 i?° 
TJ 9 > 



(2tt) c 



and 



q 4 (v k + Rl) 2 - 2 q 2 X k V^o {D k + R k ) 



d s Xk = 0. 



3/2 



(10) 

(11) 



The critical regime corresponds to a fixed point (FP) of the (dimensionless) flow of Tk, where 
the coefficients Dk and v k — denoted Xk — assume a power-law behaviour Xk ~ k~ Xx with 
x x = ~d s \n.Xk- Since, at criticality, C ~ k~ 2x and uj ~ k z , the anomalous dimensions x v 
and x D (at the FP) are related to the critical exponents x an d z through: z = 2 — x v and 
X = 1 — d/2 + {x D — x v )/2. We emphasize that d s Xk = is verified at all scales k and for all 
ipo. This implies that the Galilean invariance is preserved all along the flow at any non-zero FP, 
which in turn enforces the relation x D — 3 x v = d — 2 (or equivalently z + x = 2). Moreover, one 
can check from Eqs. © and (JTUJ) that, upon setting Dk = v k and d = 1, one obtains d s Dk = d s v k 
at all scales k, i.e. FDT is satisfied all along the flow in d = 1. 

One usually considers the combination Xk = X\Dkjv\. Upon introducing the dimensionless 
variables (p = (k d+2 i/ k /D k )~ 1 / 2 $0, 9k = k d ~ 2 X~k, y = q 2 /k 2 and the explicit cut-off function 
r(y) = (1/y — 1) 6(1 — y) [22j — where 9 denotes the Heaviside step function [2(| — we obtain: 



<9, In D k = 



-4gfe v d 



d(2 + d) m 5 / 2 



+ d + 2x D — 3x u + cp -Jgk(d + 2 — x D ^j 



d s In Vk = 



^9k Vd 



d 2 (2 + d) m 5 / 2 



d 2 - 4 - 2 x D + (4 - d) x v + (P Vfffe (rf 2 - 4 + (3 d - 2) x D - A(d - 1) x„ 
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d s gk = (d - 2) g k + g k (d s In D k - 3 d s In v k ^j (12) 

where to = 1 - 2 </> ^ and v d 1 = 2 d+1 7r d / 2 T(d/2). 

These equations are valid for arbitrary d and g k and, as we show below, they provide, at a 
qualitative level, the complete phase diagram of the KPZ equation in the (g, d) plane, which is 
displayed on Fig. [TJ This constitutes the main result of this work: the NPRG equations derived 
here entails a stable strong-coupling FP that corresponds to the rough phase in all dimensions, 
yielding in particular the exact results in d — and d = 1. 
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FIG. 1: Phase diagram of the KPZ problem. g E w, <?tr, <?p and <? K pz stands for Edwards- Wilkinson, transi- 
tion, perturbative and strong-coupling FP respectively, see text for comments. This diagram corresponds to 
4> = but remains qualitatively unchanged for non-vanishing (f>. Fhe horizontal dashed lines are guidelines 
for the eyes. 



We now present in details the FP solutions of the NPRG flow equations (fl2"|) . For cf> — 0, they 
can be solved analytically and exhibit three FP, denoted g EW , gtr and g KPZ . The FP g EW = 
corresponds to the linear Edwards- Wilkinson FP characterised by x v — x D — 0, [i.e. z = 2, 
X = 1 — d/2). The FP g tI and g KPZ are both Galilean-symmetric, i.e. satisfy z + x = 2. The FP 
gtr is unstable and hence drives a transition for d > 2. We find x v < 0, departing all the more 
from zero than the dimension increases, which indicates that the accuracy of the approximation 
deteriorates as d grows. This result means that the dynamical exponent reaches the upper 
bound z = 2 at the transition, which is consistent with the exact result x = ensuing from the 
perturbative RG re-summed to all orders [111 ]. We now focus on g K pz, which is a stable FP in 
all d. First, it yields in d = 1 the exact exponent z — 3/2. This result follows from the time 
reversal invariance ([7|) of the flow equations in d = 1. The zero-dimensional exact result z — 4/3 
is also recovered. Moreover, g K pz can be followed in all dimensions and allows to investigate 
the properties of the high-dimensional rough phase. This shows that our method is genuinely 
nonperturbative. Indeed, conversely, the perturbative one-loop RG FP g p [2l|, [23| shown on Fig. 
[T] — which also simply stems from the O(gl) expansion of Eqs. (fT2"|) — diverges in d = 3/2 where 
the perturbative scheme breaks down. Moreover the re-summed to all order perturbative RG flow 
equation for g does not yield any non-trivial stable FP, i.e. it cannot access the strong-coupling 
regime fTTj j . 

Unfortunately, while the Ansatz (|S|) is sufficient to provide the whole phase diagram at a qual- 
itative level, it does not yield correct critical exponents in d > 2: in fact the value of z starts 
decreasing for d > 2, which is unphysical. 

Let us now discuss possible improvements. A key advantage of the NPRG formalism is that it 
can be improved in a systematic way. Indeed, efficient non-perturbative approximation schemes 
have been developed and thoroughly studied during the last decade [15[. They have shown that 



the lowest order truncation (analogous to the Ansatz used here) al way s turns out to be sufficient 
to capture in a robust way all the qualitative features of a model [Tol [l7l [l8l]. For most models, 
even out of equilibrium (reaction-diffusion processes for instance) , going to the next order in the 
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derivative expansion becomes increasingly tedious but allows to obtain very accurate exponents 
(and other physical quantities) flBL [2^|. 

For the KPZ action ([5]), the next approximation level within the derivative expansion consists 
in considering field-dependent functions Vk{ip) and Dk(tp) (there is no ip dependence at this order 
due to the symmetries). We have investigated this level of truncation. Prior to discussing the 
results, let us state the outcome: unfortunately, while confirming the qualitative phase diagram 
found at the lowest order, it remains unsufficient to obtain meaningfull exponents for d > 2. 
This failure seems to root in very slow convergence properties of the derivative expansion for 
d > 2, as can be illustrated by probing the influence of the dimensionless response field value 
(j> — which is related to the noise strength. The strong-coupling FP solution of Eqs. (fT2"|) can 
be determined numerically for different finite <fi, which has been achieved. The dependence of 
gkpz on 4> m various dimensions is displayed on Fig. [2] It shows that the behaviour of g KPZ 
seems to undergo a noticeable change crossing d = 2. Below two dimensions, the FP coordinate 
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FIG. 2: Strong-coupling FP gr K pz as a function of the dimensionless response field cj>- The curves below 
(top row) and above (bottom row) d — 2 are split in two sets for clarity. Note that in d = 4, g has been 
rescaled to g/5 to fit in. 



<?kpz(</>) shows a monotonic and weak dependence on the field <f>, whereas above d — 2 it acquires 
a non-trivial shape, the two regimes appearing separated by the d — 2 curve (see Fig. [21 note the 
different scales for g). This suggests that the convergence of the scheme deteriorates above d = 2. 
In fact, the specificity of the KPZ model, compared with other non-equilibrium models such as 
reaction-diffusion processes for instance, is the nature of its interaction, which is derivative, and 
not potential. Whereas the derivative expansion converges very rapidly for models with potential 
interactions [H, [l?], HI], it would require for the KPZ model many orders to lead to quantitative 
exponents, which is not conceivable. This suggests that the derivative expansion is not the best 
strategy to investigate the physics of the rough phase. One should device an alternative scheme 
of approximation, adapted to derivative interactions. This will be presented elsewhere [25j . 

To summarize, we presented in this work the first theoretical approach to the KPZ problem, 
which entails the strong-coupling FP describing the rough phase in all dimensions, embedding 
in particular the exact results in d = and 1. While the approximation scheme used here, the 
derivative expansion, provides a consistent and robust picture of the whole phase diagram at 
a qualitative level, it is not sufficient to obtain physical critical exponents above d = 2. An 
alternative approximation scheme is needed for a thorough analysis of the properties of the 
strong-coupling regime. 
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